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We investigate in detail the model of a trophic web proposed by Amaral and Meyer [Phys. Rev. 
Lett. 82, 652 (1999)]. We focused on small-size systems that are relevant for real biological food 
webs and for which the fluctuations are playing an important role. We show, using Monte Carlo 
simulations, that such webs can be non- viable, leading to extinction of all species in small and/or 
weakly coupled systems. Estimations of the extinction times and survival chances are also given. We 
show that before the extinction the fraction of highly-connected species ( "omnivores" ) is increasing. 
Viable food webs exhibit a pyramidal structure, where the density of occupied niches is higher at 
lower trophic levels, and moreover the occupations of adjacent levels are closely correlated. We also 
demonstrate that the distribution of the lengths of food chains has an exponential character and 
changes weakly with the parameters of the model. On the contrary, the distribution of avalanche 
sizes of the extinct species depends strongly on the connectedness of the web. For rather loosely 
connected systems we recover the power-law type of behavior with the same exponent as found in 
earlier studies, while for densely-connected webs the distribution is not of a power-law type. 

PACS numbers: 87.10.+e,05.40.-a,05.45.-a,05.70.Ln 
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I. INTRODUCTION 

Food webs describe the resources and trophic relation- 
ships among species within an ecosystem. The first semi- 
quantitative descriptions of food webs were given by biol- 
ogists at the end of the nineteen century . Later on 
prey-predator relationship between species were defined 
in terms of oriented graphs with hierarchical or layered 
structures [|| . The problem of describing such food webs 
was then taken over by mathematicians and physicists, 
and different modeling levels and types of models have 
been proposed. 

A first group of models is constituted by the so-called 
static models in which the links between different species 
are assigned once and for all, according to different sce- 
narios (random, scale- free or small- world graphs 0, Q, 
for example). Some properties of these food webs were 
analyzed and compared with available biological data, 
and the comparison usually turned out to be quite poor. 

The second group of models contains the so-called dy- 
namic food web models. The novelty consists in recogniz- 
ing that the links between the species are generally not 
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arbitrary and quenched, but emerge as the result of some 
intrinsic biological dynamics. There are then many pos- 
sibilities to model the evolutionary dynamics @. The 
simplest one concerns two-layered systems with prey- 
predator Lotka-Volterra type of dynamics (for a short 
review, see Q)- A very large body of work has been 
devoted to the study of population dynamics equations 
for more than two species [1, Q. In such cases, the links 
among the species can be modified according to the evo- 
lutionary dynamics. One important issue is the control of 
the robustness of such models when the complexity of the 
system is increased. Moreover, at a more refined level of 
description, the Lotka-Volterra mean-field dynamics can 
be replaced by individual-based models [l(J, OH taking 
into account the particularities of the interacting indi- 
viduals and thus offering the possibility to include the 
stochastic fluctuations. These dynamic food webs mod- 
els allow therefore to treat on an equal footing both the 
micro- and the macro-evolution of an ecosystem fl2l [l3[ . 

The richness of the models mentioned above has its 
own drawbacks. Indeed, the number of control parame- 
ters defining the models is usually quite large; moreover 
the dynamics is nonlinear. Thus, it is often impossible 
to get a global picture of the properties of the system. 
Accordingly, it is desirable to study some models which 
are as simple as possible, in order to clarify the rela- 
tive importance of the various ingredients, while being 
able to capture the generic properties expected for food 



webs. Several proposals have been made along this line 
in the past years, see e.g. [3 [Hj]. In particular, Amaral 
and Meyer [16| proposed such a "minimal" model whose 
numerical solution leads to a power-law distribution of 
extinction-avalanche sizes, in good agreement with avail- 
able data from fossile record. It was shown later that this 
model is self-organized critical [Tt} and that the power 
law can be obtained analytically. Furthermore, taxo- 
nomic effects have been added to the model [Tj|, but 
without significant effects. 

In this work we are revisiting the Amaral - Meyer 
model (AM hereafter) with the aim of investigating sev- 
eral of its properties which are relevant for real food webs 
and which have not been addressed in the previous works. 
The paper is organized as follows. In Sec. II, the model 
is described and several technical details concerning the 
Monte-Carlo simulations, as well as the values of the con- 
trol parameters, are given. Section III contains the main 
results. First, the dependence of the survival chance and 
of the average extinction time on the number of niches TV 
and on the maximum number of feeding species k is stud- 
ied. The problem of extinction due to stochastic effects 
is also discussed. Then the question of the pyramidal 
structure of the food web is approached. Time correla- 
tions between the occupied niches at different levels are 
investigated. The time evolution of the ratio of omni- 
vores is also computed, both for viable and non-viable 
food webs. The distribution of food-tree sizes as a func- 
tion of the values of N and k is found to exhibit different 
regimes. Finally, the problem of avalanches of species 
extinctions is revisited. In contradiction with previous 
results, it is found that strong deviations from simple 
power laws for the size distribution of these avalanches 
can be observed for large values of k. Some of our pre- 
dictions are compared with real biological data and are 
found to be in good agreement. Conclusions arc relegated 
to Sec. IV. 



II. MODEL 

The AM food web model consists of L trophic levels, 
each of them containing the same number N of niches, 
which can be either empty or occupied by a single species. 
Each species from level I = 2, 3, L feeds on at most k 
(k 1) species that are randomly selected from the level 
below, (I — 1) (sec Fig.[TJ. Therefore a species from level 
I is a predator for some species at the level (I — 1), and 
at the same time it may be a prey for species from the 
level (I + 1) (except for the species on the top level L, 
that have no predators, and the species on the bottom 
level I = 1 that have no preys). 

The dynamics of the web is driven by the extinction 
and creation of species, as well as by the dynamically- 
related evolution of the trophic links between the species. 
Namely, at each time step (Monte Carlo step, MCS), 
one starts by randomly removing species from the basal 
level 1, with a given "extinction probability" p. When 
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FIG. 1: Schematic representation of the AM food web model, 
for L = 6 and N = 5. The occupied niches are represented 
by the black rectangles and the interactions between species 
are depicted by the directed links. 



a species gets extinct, all the links from it to species at 
the level 2 are removed. If as the result of these link 
removals a species at the level 2 looses all its preys from 
level 1, then it becomes extinct as well. This procedure 
of checking existing links and removal of species which 
lost all their food sources is then followed on each level 
till the top level L. Hence an avalanche of extinctions of 
species can be generated. 

Apart from extinctions, the AM model considers also 
the creation of species in the free niches. Each species 
(that remains after the decimation procedure described 
above) at level I can repopulate, with a probability /i, 
an empty niche either at the level I, (I — I), or (I + I). 
New-created species receive at most k links, at random, 
to species from the adjacent lower level. 

Extinction and creation of species are thus stochas- 
tic processes that differ from one realization of the food 
web to another, and one can address the question of 
the statistical properties of various characteristics of the 
system, like, for example, the size of the extinction 
avalanches, the extinction time (or, equivalcntly, the sur- 
vival chance), the populations at all levels, the correla- 
tions between the different trophic levels, the density of 
the trophic links, etc. The dependence of these elements 
on the parameters of the model L, N, k, p, and /i is also 
an important aspect to be considered. 

In this respect, the main result of the original pa- 
per [l(| addresing the AM model was that the distri- 
bution of the sizes of the extinction avalanches can be 
fitted over about three decades by a power law with 
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an exponent a ~ —1.98; this exponent was corrected to 
the value a — — 2 in later works [T3, EH, which is sup- 
ported by mean-field theoretical arguments. Moreover, it 
was argued that this power-law behavior is in agreement 
with available fossil data records. In Ref. it has also 
been shown that the avalanche-size distribution exhibits 
a maximum for small-size events, before dcvclopping the 
power-law behavior. However, most of the characteris- 
tics of the food web that were enumerated above were 
not addressed in the previous papers on the AM model 
and our work is therefore intended to fill this gap. 

We shall therefore investigate the AM model in more 
detail, by considering the canonical set of parameters 
used in [l6[, namely L = 7 trophic layers, the extinction 
probability p = 0.01 and the probability of creation fx = 
0.02. We shall moreover investigate how the system char- 
acteristics depend on the number of niches N and on the 
highest possible number of links k that a predator may 
have. The obtained results will be compared to experi- 
mental data coming from investigations of some contem- 
porary food webs [l9|, [2(], Hl[ . Since the total size - i.e., 
the product of the number of layers and the mean num- 
ber of occupied niches - of the experimentally-observed 
food webs does not exceed 1000, we have decided to focus 
on N values that are smaller than the value of 1000 that 
was used in Refs. [TH, [l8[, and to work with N ^ 100, 
exceptionally 200 and 500. As a consequence, the role 
of fluctuations in our systems becomes more important 
and many of the reported effects are clearly noise-induced 
and/or noise-affected, which actually makes them more 
relevant for real biological food webs. This choice of small 
N also allowed us to run the simulations for longer times 
than those considered in [HI [l8j , which unveiled new as- 
pects of the food web viability. In general, we performed 
simulations over ~ 10 6 MCS and the averaging was done 
over 100 runs (i.e., random realizations of the food- web 
stochastic dynamics). In some cases, in order to check 
the viability of the system, we even went to ~ 10 7 MCS. 
Mean extinction times for the whole web were obtained 
by averaging over 500 runs. 
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FIG. 2: Survival chance at time t for different values of N 
and k. Evaluated from 100 runs. 
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FIG. 3: Dependence of the average extinction time of the 
whoie food web on TV and k. The average was performed over 
500 runs. 



III. MAIN RESULTS 

A first result refers to the viability of the food web, 
i.e., to its capacity to survive in the long time limit. 
Performing much longer simulations than in [l8| we 
have found that small-size (e.g., N = 50 or N = 70) and 
weakly-coupled (k = 3 - 6) systems are not viable and 
disappear in the long time limit t « 10 6 MCS. Figure [2] 
illustrates how the chance that a web will survive till a 
given time t is depending on N and k. Survival chance 
at time t is defined here as the ratio of the number of 
realizations (runs) for which the system was still existing 
at time t, to the total number of trials. The web of N = 
100 and k = 3 turns out to be non- viable, too (out of 100 
runs, none has survived till 10 7 MCS), however increas- 
ing k to 4 stabilizes the system. The dependence of the 



mean extinction time on the number of niches A, and on 
the maximum number of links k is illustrated in Fig. [3] 
It is obvious that increasing k, i.e., the connectedness, is 
stabilizing the web. Small, sparsely coupled webs cannot 
exist for a longer time. 

The mechanisms leading to the collapse is connected 
with the stochastic nature of the extinction and prolifer- 
ation events. Indeed, when the system is small, it may 
happen rather easily that at the lowest level, which is 
crucial for the survival of the web, only very few species 
survive. If, moreover, as is the case in the AM model, the 
values of the two creation and extinction probabilities, 
are very low, then two scenarios are almost equally prob- 
able, namely: (i) either some empty niches at the level 1 
are repopulated and the web is, at least temporarily, safe, 
or (ii) existing species are all removed from this level, as 
illustrated in Fig. 2J This is the end of the food web, 
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FIG. 4: Temporal evolution of the number of species on the 
basal level for N = 50, k — 3. Single run, exhibiting the 
complete extinction of the species on this level. 



since without species at the basal level an avalanche con- 
taining all species is created and the web collapses. This 
stochastic extinction in small populations is a well known 
effect in ecology p2j . 

The next set of figures, see Figs. [5j illustrate the tem- 
poral evolution of the normalized populations at the dif- 
ferent levels (i.e., the number of species at a given level 
divided by the total number of species in the web, at a 
given time). For the sake of clarity, only a part of the 
levels are shown. The AM model leads in a natural way 
to a pyramidal form of the food web, where the upper lev- 
els are less populated that the lower ones, see the upper 
panel. This effect is less pronounced when the system is 
close to its collapse, as shown by the lower panel of Fig. [5] 
Comparison with Fig. [6] indicates that the pyramid-effect 
is practically disappearing for systems with many niches 
(high N), and this is the reason for which it had not 
been observed neither in [l6j|, nor in [l8[ for which N = 
1000. Note also that the pyramidal structure has been 
best observed by biologists in rather small food webs . 

As seen from both Figs. [5] and the numbers of oc- 
cupied niches at different levels are randomly oscillating 
in time. In order to check the degree of correlation of 
these oscillations, and whether there is some systematic 
time lag between them, we have calculated the correla- 
tion functions from the corresponding discrete-time series 
of the populations on the different levels, using the for- 
mula [23|: 

( T-m 

^2 SN i ( n + m ) sn j ( n ) > for m > 0, 



Cij(m) = C t < 



n=0 



T-\m\ 



^2 SNi(n) SNj(n + \m\) , for m < 

I n=0 

Here iVj (t) is the population of level I = i at time t (which 
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FIG. 5: Time dependence of the normalized populations of 
different levels for N = 100, k = 3. Upper panel: initial and 
mean time stage. Lower panel: long-time stage, just before 
the web collapse. Single runs. 



is, of course, an integer number of MCS) and SNi(t) is 
its fluctuation around the mean value, 



5Ni(t) = Ni(t) - T- 1 SNi(n); 

n=0 

m is the time lag (that can be positive or negative), and 
T is the total simulation time. The coefficient C,; was 
chosen such that the autocorrelation functions at zero 
lag are equal to 1, 



The results for N = 100 and N = 200 are illustrated in 
Figs. [7] and [HJ respectively. From these figures it fol- 
lows that the time series at neighboring levels are highly 
correlated with each other at zero time lag. In other 
words, species at a given level adjust immediately to the 
changes at the level below, which is a feature that could 
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FIG. 6: Time dependence of the normalized populations of 
different levels for k = 3. Upper panel: TV = 200. Lower 
panel: TV = 500, single runs. 
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FIG. 7: Cross-correlations between the occupation numbers 
at levels i = 1 and j = 2, 4, 6. Here TV = 100 and k = 3, and 
the evaluations are done on a single time-series. 
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FIG. 8: Cross-correlations between the occupation numbers 
at levels i = 1 and j = 2, 4, 6. TV = 200, fc = 3, evaluations 
from single runs. 



be expected in view of the constitutive dynamics of the 
model. The correlation is of course decreasing with the 
distance between the levels, but the peak at zero time 
lag remains. It should be noted that in the TV = 100 
case, Fig. [JJ apart from the very narrow zero-lag correla- 
tion peak there is also a rather broad structure centered 
around it. This structure is practically absent in the case 
of the larger system with TV = 200, sec Fig. [5] and is to 
be related to the long-time instability of the system with 
TV = 100 and k = 3, and to the strong fluctuations that 
are accompanying its collapse. 

Another biologically interesting feature is the frac- 
tion of omnivores, that are predators feeding on more 
than one prey [24| ■ Figures [9] shows the distribution of 
the number of links per predator for a non-viable (up- 
per panel), and a viable (lower panel) system, at sev- 
eral times. The distribution remains virtually the same 
throughout the simulation time for a viable web. How- 
ever, for a non-viable one approaching extinction, the 
fraction of highly-connected predators grows. In other 
words, close to the collapse, only predators feeding on 
many preys will survive. In this sense, the presence of 
omnivores stabilizes the web, as documented experimen- 
tally in [IJ. 

The average fraction of omnivores in a stationary state 
of the food web depends on the maximum value of links 
k, but it is rather insensitive to the number of niches TV, 
as illustrated by the data in Table HI This feature also 
agrees with experimental results reported in [2lj . 

In a stationary state of the web, the average number 
of links corresponding to different levels I has the same 
type of profile whatever the value of k, namely a more 
or less pronounced maximum for the intermediary levels 
and a drop for the low and top levels, as illustrated in 

Fig.nni 

Food webs are also often characterized by the lengt h o f 
the food chains (or "trees" ) that are forming the web [19l . 
[2l| . We define them in the following way: each species 
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FIG. 9: Distribution of the number of links per predator for 
non- viable k — 6 (upper plot), and viable k = 20 (lower plot) 
systems, at different times. N = 50 in both plots, evaluations 
from single runs. 
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FIG. 10: Number of links for different levels in a stationary 
state of the web, averaged over 100 runs. 
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FIG. 11: Food tree size distribution at time T = 3 • 10 5 MCS 
as a function of the number of niches TV. k = 6, semi- 
logarithmic scale, evaluation from 100 runs. 



TABLE I: Average fraction of omnivores in a stationary state 
of the web, for various k and TV. The average was taken over 
100 runs. 



with no predators is the root of a new tree. Starting 
from the root we go along its links to the lower level and 
mark all species the root is predating on. Then we check 
their links to find their prey species and so on. Since 
different predators do not really compete for food in the 
model (i.e., if they are linked to the same prey, they all 
get enough food), we can treat the partially overlapping 
trees as independent ones. The size of a tree is then 
simply the total number of species that belong to that 
tree. 

As can be seen from Fig. [TTJ food tree size distribution 



depends on the number of niches N in the system. The 
maximal tree size increases with N, as could be expected. 
Moreover, in a bigger system there is more space for trees 
of similar sizes and that is why the curves in Fig. [TT] shift 
upwards with increasing N. Linear dependence of the 
distribution of chain lengths on the semi-logarithmic plot 
in Fig. [TTJ indicates an exponential decrease with the tree 
size. 

The maximum number of links k between species is 
also playing an important role on the food tree size dis- 
tribution. The results for AT = 100 niches at two different 
time steps are displayed in Fig. [T2] When k increases, 
small trees become less likely and bigger structures in 
the system are preferred instead. We can thus distin- 
guish two regimes with different ^-dependence, namely 
the regime of "small trees" (of size ^ 10), whose num- 
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FIG. 12: Food tree size distribution at time T — 5 ■ 10 4 MCS 
(upper plot) and T = 3 • 10 J MCS (lower plot) as a function 
of the maximal number of links between species k. N — 100, 
log-log scale, evaluation from 100 runs. 
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FIG. 13: Distribution of avalanche sizes for N = 200 and N 
= 500, for food webs with k = 3. Estimated from 500 runs. 
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FIG. 14: Distribution of avalanches for N = 100 and different 
k. Estimated from 500 runs. 



TABLE II: Fraction of species at basal, intermediary and top 
levels, k — 6, average over 100 runs. 



ber decreases with k, while the number of "big trees" 
is an increasing function of k. Moreover, the maximal 
size of a food tree varies strongly with k. It is also in- 
teresting to note that for sufficiently large N there is a 
well-pronounced peak in the distribution of tree sizes at 
a size of 7, which is simply the number of trophic levels 
in the system. 

Biologists [l], [2l[ often describe food webs in terms of 
fractions of basal, intermediate and top species. In this 
model these ones correspond to occupied niches at levels 
1, 2-6, and 7, respectively. For viable systems the values 
we obtain for these fractions are not too sensitive to the 
values of N and are presented in Tabic [Til These results 



agree very well with biological data for a food web from 
Little Rock, U.S.A., see [2l| for further details. 

Finally, we have analyzed the distribution of avalanche 
sizes of species extinctions. We have observed the max- 
imum in the distribution that was mentioned in [l8l ]. 
which becomes more pronounced with increasing the 
number of niches. For k — 3, which was the value 
considered in [lj| [l8| , we recovered the known power-law 
behavior, extending over nearly three decades, with an 
exponent equal to a — —2, as calculated in [Tt], EU- This 
value does not seem to depend on N, and even for non- 
viable systems we got the same good fit to a power-law 
type of behavior, with the same exponent, see Figs. 1131 

EE3 

However, when the food web becomes highly connected 
(i.e., k = 6 or larger) the deviations from the power-law 
behavior are very large, as illustrated in Figs. RBI IT51 for 
k = 20. One may notice that for highly-connected webs 
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the fraction of larger avalanches increases, simply because 
the removal of a prey on which many predators feed is 
affecting more species. 

As such, the power-law behavior, although largely 
present in the large-A webs, is not universal, but it is 
rather the result of a particular choice of the parameter 
values of the food web. 



IV. CONCLUSIONS 

We have presented a detailed discussion of several as- 
pects of the AM model [n| of a food web that were ne- 
glected in the previous studies. In particular, we concen- 
trated on relatively small-size systems and on the role 
the fluctuations can play in such systems, since this is 
the frame that is important in most of the real biolog- 
ical food webs. Several comparisons of the theoretical 
predictions with experimental data were also discussed. 

Our simulations confirmed the observations of [l8j con- 
cerning the distribution of avalanche sizes of species ex- 
tinctions, the value of the power-law exponent a = — 2, 
and the existence of a maximum depending on N. It is 
worth noting here that the same value a = — 2 of the 
exponent of the distribution of extinction avalanches has 
also been found in two other different food web models, 



Refs. [25| and [26|. The last paper is a generalization of 
the Bak and Sneppen model [271 ]. in which two factors 
determine the fate of a species - biotic ( "bad genes" ) and 
abiotic ("bad luck"). 

For reference, we have kept the values of some of the 
model parameters (number of layers L, probability of cre- 
ation \x and of extinction p) the same as in the original 
AM model [lj| [H| . However the dependence of the sys- 
tem behavior on the number of niches, N and on the 
maximum number of links per predator k, turned out to 
be quite interesting. We have thus unveiled new features 
of the model, not found in the earlier papers. A food web 
may collapse if it is too small and/or has not enough 
links between species. Systems smaller than N 200 
show a pyramid-like structure, where top levels are less 
populated than the bottom ones. The occupations of the 
levels are strongly correlated at zero time lag. When the 
web is close to a collapse, the fraction of highly connected 
predators (omnivores) significantly increases, which may 
lead sometimes to a (temporary) rescue of the web. The 
distribution of the length of food trees has an exponen- 
tial character and its type depends rather weakly on both 
N and k. Finally, the distribution of species extinctions 
shows an unexpected feature, contrary to the previously- 
claimed universality of the power-law behavior; namely 
that for large A:- values (i.e., highly-connected webs) the 
distribution cannot be described anymore as a power law. 

Although several criticisms concerning the applicabil- 
ity of the AM model to biology have been raised [l8l . |28| , 
we have found that some theoretical results, like the ratio 
of omnivores, the fraction of different-type (basal, mid- 
dle, top) species, the food chain length etc., are in very 
good agreement with experimental data [2(| Hl[ on 
food webs. The significance of this fact is a subject of 
further analysis. 



Acknowledgments 

This work has been done within the UNESCO Chair 
for Interdisciplinary Studies activity. AP and JS ac- 
knowledge financial support from COST P10 Project 
"Physics of Risk" and the hospitality of the Departcment 
de Physique Thcorique de l'Universite de Geneve. This 
work was partly supported by the Swiss National Science 
Foundation. 



[1] R. T. Paine, J. Anim. Ecol. 49, 667 (1980). 

[2] J. E. Cohen, F. Briand, and C. M. Newman, Biomathe- 

matics Vol. 20, Springer- Verlag Berlin, 1990. 
[3] S. L. Pimm, Food Webs, Chapman & Hall, London, 1982. 
[4] S. H. Strogatz, Nature 410, 268 (2001). 
[5] R. Albert and A. L. Barabasi, Rev. Mod. Phys. 74, 47 

(2002). 

[6] B. Drossel and A. McKane, in Handbook of graphs and 



networks, S. Bornholdt and H. G. Schuster (Edts.), 

Wiley- VCH, Berlin, 2003. 
[7] A. Pekalski, Comp. Sci. Eng. 6, 62 (2004). 
[8] R. M. May, Stability and complexity in model ecosystems, 

Princeton University Press, 1974. 
[9] B. Drossel, Adv. Phys. 50, 209 (2001). 
[10] C. Quince, P. G. Higgs, and A. J. McKane, Ecol. Model. 

187, 389 (2005). 



9 



[11] M. Droz, Individually based models approach to the dy- 
namics of cooperative behavior in social sciences, in Mod- 
eling Cooperative Behavior in the Social Sciences, P. L. 
Garrido, J. Marro and M. A. Munoz (Edts.), AIP con- 
ference proceedings 779, 29 (2005). 

[12] D. Chowdhury and D. Stauffer, Phys. Rev. E 68, 041901 
(2003); D. Stauffer and D. Chowdhury, Physica A 336, 
102 (2004). 

[13] F. Coppex, M. Droz, and A. Lipowski, Phys. Rev. E 69, 

061901 (2004). 
[14] G. Abramson, Phys. Rev. E 55, 785 (1997). 
[15] R. V. Sole and S. C. Manrubia, Phys.Rev. E 54, R42 

(1996); ibid, Phys. Rev. E 55, 4500 (1997). 
[16] L. A. N. Amaral and M. Meyer, Phys. Rev. Lett. 82, 652 

(1999). 

[17] B. Drossel, Phys. Rev. Lett. 81, 5011 (1998). 

[18] J. Camacho and R. V. Sole, Phys. Rev. E 62, 1119 (2000). 

[19] C. R. Townsend et al, Ecol. Lett. 1, 200 (1998). 



[20] J. A. Dunne, R. J. Williams, and N. D. Martinez, Ecol. 

Lett. 5, 558 (2002). 
[21] R. J. Williams and N. D. Martinez, Nature 404, 180 

(2000). 

[22] M. Shaffer, BioScience, 31, 131 (1981). 

[23] S. J. Orfanidis, Optimum Signal Processing. An Intro- 
duction., 2nd Edition, Prentice-Hall, Englewood Cliffs, 
1996. 

[24] C. Borrval, B. Ebenman, and T. Jonson, Ecol. Lett. 3, 
131 (2000). 

[25] P. A. Rikvold and R. K. P. Zia, Phys. Rev. E 68, 031913 
(2003). 

[26] B. W. Roberts and M. E. J. Newman, J. Theor. Biol. 
180, 39 (1996). 

[27] P. Bak and K. Sneppen, Phys. Rev. Lett. 71, 4083 (1993). 
[28] M. E. J. Newman and R. B. Palmer, Modeling extinction, 
Oxford University Press, Oxford, 2002. 



